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ABSTRACT 

Context. The space experiment CoRoT has recently detected transits by a hot Jupiter across the disk of an active G7V star (CoRoT- 
Exo-2a) that can be considered as a good proxy for the Sun at an age of approximately 0.5 Gyr 

Aims. We present a spot modelling of the optical variability of the star during 142 days of uninterrupted observations performed by 
CoRoT with unprecedented photometric precision. 

Methods. We apply spot modelling approaches previously tested in the case of the Sun by modelling total solar irradiance variations, 
a good proxy for the optical flux variations of the Sun as a star. The best results in terms of mapping of the surface brightness 
inhomogeneities are obtained by means of maximum entropy regularized models. To model the light curve of CoRoT-Exo-2a, we 
take into account both the photometric effects of cool spots as well as those of solar-like faculae, adopting solar analogy. 
Results. Two active longitudes initially on opposite hemispheres are found on the photosphere of CoRoT-Exo-2a with a rotation 
period of 4.522 ± 0.024 days. Their separation changes by x 80° during the time span of the observations. From this variation, a 
relative amplitude of the surface differential rotation lower than ~ 1 percent is estimated. Individual spots form within the active 
longitudes and show an angular velocity about ~ 1 percent smaller than that of the longitude pattern. The total spotted area shows a 
cyclic oscillation with a period of 28.9 ± 4.3 days, which is close to 10 times the synodic period of the planet as seen by the rotating 
active longitudes. We discuss the effects of solar-like faculae on our models, finding indication of a facular contribution to the optical 
flux variations of CoRoT-Exo-2a being significantly smaller than in the present Sun. 

Conclusions. The implications of such results for the internal rotation of CoRoT-Exo-2a are discussed on the basis of solar analogy. 
A possible magnetic star-planet interaction is suggested by the cyclic variation of the spotted area. Alternatively, the 28.9-d cycle may 
be related to Rossby-type waves propagating in the subphotospheric layers of the star. 

Key words, stars: magnetic fields ~ stars: late-type ~ stars: activity - stars: rotation - planetary systems - stars: individual (CoRoT- 
Exo-2a) 



1. Introduction 

The CoRoT satellite is a space experiment devoted to asteroseis- 
mology and extrasolar planet search by the observa tion of plan- 
etary transits across the discs of their parent stars CBaglin et alJ 
|2006|) . One of the first planets detected is CoRoT-Exo-2b, a hot 
Jupiter orbiting a main-sequence G7 star with a period of 1 .743 



* Based on observations obtained with CoRoT, a space project op- 
erated by the French Space Agency, CNES, with partecipation of the 
Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, 
Germany, and Spain. 



days. The planet's characteristics and the physical parameters 
of the star have been determined by com bining CoRoT pho- 
tometry with groun d-based spectroscopy jAlonso et al.l 120081 : 
iBouchv et al.ll2008h . In the present work we address the mag- 
netic activity of the star as revealed by the rotational modulation 
of its broad-band optical flux induced by photospheric bright- 
ness inhomogeneities. The short rotation period (~ 4.5 days) 
and the sizeable convective envelope of CoRoT-Exo-2a promote 
a vigorous dynamo action that amplifies and modulates its mag- 
netic fields. Magnetic flux tubes emerging to the surface give 
rise to cool spots and bright faculae in the photosphere of this 
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Star whose mass, radiu s and effective temp erature are compara- 
ble to those of the Sun. lBouchv et alj ( l2008l) suggest that the star 
is still close to the ZAMS, possibly younger than 0.5 Gyr, repre- 
senting a unique laboratory to study the magnetic activity of the 
young Sun, at a level of flux variations one order of magnitude 
greater than the present Sun. 

We present the results of a spot modelling of CoRoT-Exo-2a 
based on the methods already applied to model the bolometri c 
variabihty of the Sun as a star jLanza et alJl2003L |2004 l2007l) . 
The uninterrupted time series and the high photometric precision 
of CoRoT observations allow us to map the longitudinal distri- 
bution of the active regions in the photosphere of the star and fol- 
low their evolution continuously along a period of five months 
which has not been previously possible with ground-based ob- 
servations. Our approach is based on the modelling of stellar 
flux variations outside transits. The possibility to map the spots 
located on the band of the star occulted by the planet during 
its transits, by means o f eclipse mapping techniques (e.g. ISilval 
120031: iPont et al.ll2007h . wiU be the subject of another investiga- 
tion that will benefit from the present study to put its results in a 
broader context. 



2. Observations 

CoRoT-Exo-2a was observed during the first long run from May 
16 to October 5, 2007 with a time sampling of 512 s during the 
first week and 32 s thereinafter CoRoT is on a low-E arth po- 
lar orbit with a period of 1.7 hr (lAuvergne et al.ll2008l) . so the 
observations during the passage of the satellite across the South 
Atlantic Anomaly of the Earth magnetosphere are systematically 
afifected by proton hits and must be discarded. These inteiTup- 
tions usually do not exceed 15-20 percent of each orbit and 
have no significant impact for our study since the modification 
of the photospheric pattern of an active star occurs on timescales 
of the order of one day or longer. Data were corrected for the 
CCD zero oflrset and gain. The contribution of the background 
light was estimated and coiTected by observing a window of 100 
pixels close to the target. After reduction and correction for the 
known instrumental effects, those observations provide us with 
one of the highest quality CoRoT light curves. 

The passband of the photometric data used in the present 
study ranges from 300 to 11 00 nm. We preferred to use CoRoT 
white light data instead of the three colours because of the 
significantly higher signal-to-noise ratio, reduced impact of 
the hot pixels, and long-term stability of the signal. Transits 
were removed by means of the ephemeris and parameters of 
lAlonso et al.l (2008). To eliminate residual outliers, we disre- 
garded all data points that deviated from the mean more than 
4.2 standard deviations of the whole data set. Such a broad limit 
was adopted to initially include data points that have a deviation 
from the mean comparable with the depth of the transits. Then 
we computed a filtered version of the light curve by means of 
a sliding median boxcar filter with a boxcar ex tension equal to 
one o rbital period of the satellite, i.e., 6184 s (cf. lAuvergne et alj 
12008 ). This filtered light curve was subtracted from the original 
light curve and all the points deviating more than 3 standard de- 
viations of the residuals were discarded. In this way, 0.77 per- 
cent of the data initially validated by the CoRoT N2 reduction 
pipeline were discarded. Finally, we computed normal points by 
binning the data on time intervals having the duration of the or- 
bital period of the satellite, getting a light curve consisting of 
1945 normal points that cover 142.006 days. Each normal point 
is obtained by averaging ~ 193 observations with a 32-s sam- 



pling and its mean standard error is 1.03 x 10"^ in relative flux 
units. 

As pointed out bv lAlonso et al.l ( l2008l) . a faint apparent com- 

panion located 4 apart falls completely inside the photometric 
mask of CoRoT-Exo-2a contributing to 5.6±0.3 percent of the 
total flux. This fraction of the median flux has been subtracted 
and the light curve normalized to its maximum value of 701846 
(e"/32 s), observed at HJD = 2454374.2861, that we assume to 
represent the unspotted flux level of the star, whose true value is 
unknown. 

3. Spot modelling 

To fit the high precision light curves of solar-like stars provided 
by space-borne photometers, we apply two approaches that have 
been tested using obse rvations of the tota l solar irradiance, i.e., 
of the Sun as a star, bv lLanza et alj (|2007[) . The first is based on 
the use of three active regions, containing both cool spots and 
warm faculae, to fit the rotational modulation of the flux, and of 
a uniform component to fit the variation of the mean light level 
produced by a uniform background of active regions. This model 
will be hereinafter referred to as the 3 -spo t model and it is de- 
scribed in detail in Lanza et alj ([2003, 2007). Its main advantage 
is the limited number of free parameters. When the inclination 
of the stellar rotation axis is fixed, as in our case (see Sect.|4]l, we 
fit the model to the observed flux variations by varying the area, 
longitude and latitude of the three active regions, their rotation 
period and the mean flux level, having in total 1 1 free param- 
eters. The limb darkening parameters and the spot and facular 
contrasts are dete rmined by fixing the stel lar, spot and facular 
temperatures (see Lanza et al.ll2004 l2007l) . The ratio Q of the 
facular to the spotted area in the active regions is determined by 
minimizing the of the best fit t o the en tire light curve, accord- 
ing to the method suggested bv .GondoinI ( l2008i) (see Sect.|4|. 

A second method applied to fit the hght curves of solar-like 
stars is based on the use of a continuous distribution of active re- 
gions to reproduce the complex surface patterns observed on the 
Sun close to the maximum of the eleven-year cyc le or on stars 
that a re more active than the Sun itself (cf., e.g., iKovari et al.l 
|2004|) . The star is subdivided into a large number of surface el- 
ements, in our case 200 squares of side 18°, with each element 
containing unperturbed photosphere, cool spots and facular ar- 
eas. The fraction of an element area covered by cool spots is 
indicated by the filling factor /, the fractional area of the faculae 
is Qf and the fractional area of the unperturbed photosphere is 
1 - (Q + 1)/. We fit the light curve by changing the value of / 
over the surface of the star, while Q is held uniform and constant. 
Even fixing the rotati on period, the inch nation, and the spot and 
facular contrasts (see ^ Lanza et"ani200 7', for details), the model 
has 200 free parameters and sufl'ers from non-uniqueness and 
instability. To find a unique and stable spot map, we apply max- 
imum entropy reg ularization (hereinafter ME), as described in 
iLanza et a n (l2007h . by minimizing a functional which is a lin- 
ear combination of flie;if^ and the entropy functional 5, i.e.: 

@^x\f)-'^S(f), (1) 

where / is the vector of the filling factors of the surface ele- 
ments, A > a Lagrangian multiplier determining the trade-off" 
between light c urve fitting and regu larization, and the expression 
of S is given in lLanza et al.l (Il998 b. To fix the optimal value of 
the Lagrangian multiplier, we compute the mean of the residual 
for the regularized best fit yUieg and compare it with the standard 
error obtained for the unregularized best fit, i.e., eo - ctq/ y/N, 
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where cry is the standard deviation of the residuals obtained for 
A = Q and is the number of normal points in the data set. We 
fix the optimal value of A by requiring that Ifiregl - eo, i-C-, we 
iterate on A until the mean of the residuals of the ME best fit is 
equal in modulus to one standard error of the residuals of the 
unregularized model. 

Given the evolution of the sunspot configuration on the Sun, 
it is possible to obtain a good fit of the irradiance changes, both 
with the 3-spot and the ME models, only for a limited time inter- 
val Atf, not exceeding 14 days that is the lifetime of the largest 
sunspot groups dominating the irradiance variation. In the case 
of other active stars, the value of Atf must be determined from 
the observations themselves, looking for the maximum data ex- 
tension that allows for a good fit with the applied models (see 
Sect.lHfor CoRoT-Exo-2a). 

The activity level of CoRoT-Exo-2a is significantly higher 
than that of the Sun, as indicated by the amplitude of the flux 
modulation which is ~ 20 times greater than th at of the Sun 
at the maximum of the eleven-year cycle (cf. lAlonso et al.l 
120081) . Therefore, we expect that its light variations are domi- 
nated by cool spots, as suggested by Radick et al. ( 1998) and 
iLockwood et al.l (l2007h . In fact, the light curve variations can be 
accurately modelled by considering only cool spots, as we shall 
show in Sect.|5] but the possible role of solar-like faculae cannot 
be excluded a priori, so we shall consider also models with a fac- 
ular component. In Sect. |4] we shall use the 3-spot model, that 
has a number of free parameters much smaller than that of the 
ME model, to estimate the best value of the facular-to-spotted 
area ratio Q. 

The inclination of the rotation axis of CoRoT-Exo-2a is al- 
most 90° (Alonso et al. 2008), thus the duration of the transit of 
a given spot across the stellar disc is independent of its latitude in 
the hypothesis that surface differential rotation is negligible. In 
the case of the Sun, the angular velocity diff'erence produced by 
the latitudinal differential rotation in the sunspot belt (i.e. ±40° 
from the equator) is only ~ 4 percent, i.e. too small to be ap- 
preciated given the limited lifetimes of the sunspot themselves. 
Therefore, even if the amplitude of the differential rotation in 
CoRoT-Exo-2a were comparable to that of the Sun, we do not 
expect an observable dependence of the spot transit time on the 
latitude. In other words, no information on spot latitude can be 
extracted from the light curve modelling. On the other hand, the 
spot longitudes and the variation of the total spotted area during 
the 142 days of observations can be reliably extracted from the 
light curve. Our extended tests conducted on the total solar irra- 
diance variation along cycle 23 show that ME models are supe- 
rior to the 3-spot models to derive the distribution of the spotted 
area vs. longitude and the variation of the total spotted area. The 
longitude of the largest spot groups can be determined with an 
accuracy of at least 15° - 20° and the resolution in longitude, 
i.e., the capability of separating neighbour groups of compara- 
ble ai-ea, is of ~ 50° - 60° (cf. iLanza et alj|2007h . If solar-like 
faculae contribute to the light modulation of the star, the centre- 
to-limb variation of their brightness contrast may systematically 
aff bct the derivation of the active region longitudes, as discussed 
by iLanza et al. I (I2007h . Therefore, another motivation for com- 
puting models with a facular contribution is to estimate the sys- 
tematic errors of active region longitudes. 

4. Model parameters 

The basic stellar parame ters are taken from lAlonso et al.l (l2008l) 
and lBouchv et all ( l2008l) and are liste d in Tabl jll We adop t the 
same quadratic limb-darkening law of lAlonso et al.l (l2008l) with 
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Fig. 1. The ratio of the of the composite best fit to the en- 
tire time series of CoRoT-Exo-2a to its minimum value vs. the 
parameter Q, i.e., the ratio of the area of the faculae to that of 
the cool spots in active regions. The horizontal dashed line in- 
dicates the 99 percent confidence level fox lx\{m^ determining 
the interval of the acceptable Q values. 

the parameters derived from the fitting of the planetary transit 
profile. The inclination of the stellar rot ation axis is constrai ned 
by the Rossiter-McLaughlin eff'ect (see iBouchy et alJl2008l) . so 
we can assume that it is perpendicular to the orbital plane of 
the planet. We assume an isophotal wavelength of the CoRoT 
passband of 500 nm to take into account the surface distortion 
and the gravity darkening effect induced by the centrifugal force. 
Considering a rotation period of 4.52 days (see below) the rela- 
tive difference between the polar and the equatorial axes e, eval- 
uated with the Roche model, is 2.48 x 10 which may induce 
a flux variation about one order of magnitude smaller for a spot 
coverage of 10 percent of the stellar surface. The time-dependent 
tidal distortion due to the close-by planet, having a mass of 3.3 
Jupiter masses, has been computed with the model of lPfahl et al.l 
(2008) and induces a negligible flux variation with a relative am- 
plitude of ~ 3 X 10"^. The flux variation associated with the 
Doppler shift produced by the star' s reflex motion has a rela- 
tive amplitude of only 8.0 X 10"^ jLoeb & Gaudill2003h . The 
light reflected by the planet reaches 1.0 x 10 in relative flux 
units if we adopt a geometric albedo of 13 percent, that is quite 
high in view of the hot Jupi t er atmospheri c mode l s computed by, 
e.g.. ISudarskv et all ( l2000h . ISeager et all ( l2000h . iBurrows et all 
(I2008ir They are supported by t he upper Uni i t for the alb edo of 
HD 209458b recently derived bv'Rowe et al.l (l2006ll2008l) . 

The rotation period has been found by mean s of the Lomb- 
Scargle periodogram (lHorne&Baliunasl ll986') and IS -Trot — 
4.52 + 0.14 days, with the uncertainty dominated by the finite 
time extension of the data set. Since starspots are the tracers used 
to find the rotation period and considering that they evolve and 
migrate in longitude versus time, a refinement of the value given 
by the Lomb-Scargle periodogram can be found by means of the 
spot modelling itself. After some tests, we found that a value 
of the rotation period of 4.5221 days minimized the longitudi- 
nal migration of the active longitudes where new starspots were 
forming during the interval of the observations, thus we adopted 
this value as the rotation period for our modelling (cf. Sect.jS}. 

The maximum time interval Aff that can be fitted both with 
the 3-spot and the ME models turned out to be ~ 3.2 days. We 
adopted a value Aff = 3.15611 days in order to have 45 equal 
intervals covering the entire time series. Note that even if each 
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Fig. 2. Upper panel: Out-of-transit observations of CoRoT-Exo-2a (filled dots) versus time with the best fit obtained with the ME 
model including only cool spots (solid line); lower panel: the residuals of the best fit in relative flux units versus time. 



fitted data subset is only ~ 70 percent of a rotation period, the re- 
sults obtained with the ME spot modelling are fully meaningful, 
as confirmed by the tests performed in the case of the Sun where 
the models were computed using subset s with an extensio n of 
only ~ 50 percent of a solar rotation (cf. iLanza et"ani2007l for 
details). 

The brightness contrast of the spots - 0.665 was as- 
sumed to be the sam e of the sunspot bolometric contrast 
dChapman et al.l IT994h and i s independen t of the limb angle. 



as discussed by ' Lanza etaTI (I2003L |2004|) . It corresponds to a 
spot temperature ~ 540 K below that of the unperturbed pho- 
tosphere that, in the case of the Sun, accounts for the effects of 
the extended penumbra in large sunspot groups. Considering 
stars with photom etric variations of a few 0.1 magnitudes, 
iBerdvuginal ( |2005 |) found a relationship between the temperature 
of the spots and that of the unperturbed photosphere. In the case 
of CoRoT-Exo-2a, it would imply a temperature deficit as large 
as ~ 1700 K, corresponding to Cj ~ 0.21. The real temperature 
deficit of the spots on our star is probably somewhere in between 
the extreme values considered above, given that its level of ac- 
tivity is intermediate between that of the Sun and those of the 
most active stars. We prefer to adopt a strict solar analogy for 
the present analysis assuming that the active regions of CoRoT- 



Exo-2a consist of spots analogous to sunspots, but with a larger 
area or filling factor. However, we briefly discuss at the end of 
Sect.|5]how our results are affected if the spot contrast is changed 
or is not uniform over the photosphere of the star. 

In the case of the models including the facular contribution, 
we assume the same facular contrast of the Sun Cf - 0.115, with 
the same dependence on the limb angle, giving a zero contrast 
at the centre of th e disc and a maximum contrast at the limb 
(ILanza et al.ll2004l) . 

The best value of the area ratio Q between the faculae and 
the spots in each active region has been estimated by means of 
the 3-spot model. In Fig.[Tl we plot the mXio I xlnm '-^^ ^o^^\ 
of the composite best fit to the entire time series to its mini- 
mum value , versus Q. The horizontal dashed line indicates 

^ mm' *^ 

the 99 percent confidence level as derived from the F-statistics 
(e.g., Lamp ton et al. 1976). The best value of Q turns out to 
h& Q — 1.5, with an acceptable range extending from ~ 1 to 
~ 2.5. For comparison, the best value in the case of the Sun is 
Qq - 9, indicating a lower relative contribution of the faculae 
to the total light variation in late-type stars more active than the 
Sun, as suggested by Gondoin (2008). Considering the colour 
light curves of CoRoT-Exo-2a, we searched for the character- 
istic double light maxima produced by the transits of faculae 
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Table 1. Parameters adopted for the spot modelling 



Parameter 




Ref." 


Star Mass (Mq) 


0.97 


A08 


Star Radius (Rq) 


0.90 


A08 


r.ff (K) 


5625 


B08 


Limb darkening 


0.41 


A08 


Limb darlcening u^, 


0.06 


A08 


Prot (d) 


4.5221 


L08 


e 


2.48 X 10-'' 


L08 


Inclination (deg) 


87.84 


A08 


Cs 


0.665 


L08 


Cf 


0.115 


L04 


Q 


0.0, 1.5 


L08 


Aff (d) 


3.15611 


L08 



" References: A08: Alonso et al. (2008); BOS: Bouchy et al. (2008); 
L04: Lanza et al. (2004); L08: present study. 
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Fig. 3. Distribution of the residuals of the best fits of the CoRoT- 
Exo-2a light curve; dashed line: model with cool spots and with- 
out regularization (/I = 0); solid line: ME regularized model with 
only cool spots; dash-dotted line: ME regularized model with 
cool spots and solar-like faculae. The vertical dotted Une marks 
the zero value of the residuals. 

across the disc of the star, particularly evident at shorter wave- 
lengths, as it is often observed in the Sun. We found no evidence 
of such facular signatures, confirming that the facular contribu- 
tion, if present, is indeed significantly smaller than in the case of 
the Sun. This may suggest that for a low filling factor of the mag- 
netic flux tubes, facular brightness enhancement would domi- 
nate, whereas when the filling factor exceeds a certain threshold, 
convective heat transport is inhibited and dark spots appe ar with 
an inc reasingly relative contribution (cf. also iSolanki & Unruhl 
12004 . 

5. Results 

The sequence of best fits obtained with the ME model including 
only the photometric eff'ects of spots is shown in Fig.|2] together 
with the residuals. 

The distribution of the residuals is shown in Fig. |3] for the 
cases with cool spots without any regularization (A = 0), with 
only cool spots and ME regularization, and with spots and fac- 
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ulae and ME regularization. The model with only cool spots 
and without regularization has a mean of the residuals /ijes = 
1 .27 X lO"*" and a standard deviation ctq - 1 .44 x 10"'* in relative 
flux units. The model with regularization and only cool spots 
has jUies - -3.29 x 10"^ and cr - 2.26 x lO"**, that including also 
solar-like faculae has fi^es - -3.32 x 10"^ and cr = 2.39 x 10"^ in 
relative flux units. Note that the values of the Lagrangian mul- 
tipliers have been fixed by the condition \fires\ = o"o/ where 
~ 45 is the number of observations in each subset of duration 
Aff = 3.15611 days. 

The distribution of the spot covering factor / versus longi- 
tude and time is plotted in Fig. |4] for the ME models with spots 
only, and in Fig. |5] for the ME models with spots and faculae, 
respectively. In both cases the longitude increases in the same 
direction of stellar rotation and the orbital motion of the planet. 
In the case of the model with spots only, we clearly see that spots 
form within two active longitudes initially separated by ~ 180°. 
The one centred around longitude 0° does not appreciably mi- 
grate during the ~ 140 days of the observations showing a rota- 
tion period of 4.522 + 0.024 days, i.e., the period adopted for our 
spot modelling, while the other shows a retrograde migration of 
about a! 80° during the same interval, which corresponds to a 
rotation period of 4.554 days. If such a difference corresponds 
to the angular velocity at different latitudes on a differentially 
rotating star, the relative amplitude of the diff'erential rotation 
would be only 0.7 percent. Note that the angular velocity dif- 
fe rence between pole and equator derived from the fit in Fig. 2 
of iBarnes et al.l (|2005|) is about one order of magnitude greater, 
suggesting that the spots observed on CoRoT-Exo-2a are local- 
ized in a quite naiTow latitude band or that our star has been 
observed in a state with a low degree of differential rotation (cf. 
iDonati et al.l l2003; Lanza 2006i). 

The formation and evolution of individual active regions 
within each of the active longitudes is clearer in Fig.|4] It shows 
that an active region forms around HJD 2454245 in the ac- 
tive longitude centred at 0° and then increases its area and mi- 
grates backward in longitude, reaching a maximum area at HJD 
2454270, and eventuaUy fading away at HJD 2454300, after 
~ 55 days since its appearance. At HJD 2454275, a new ac- 
tive region forms in the other active longitude and evolves in a 
similar way, migrating backward in longitude and fading away 
at HJD 2454330, again with a characteristic lifetime of about 55 
days. Between HJD 2454330 and HJD 2454350, both active lon- 
gitudes show spots with similar area that decay in 20 - 30 days, 
leaving a minimum spotted area at HJD 2454350, when a large 
active region begins to develop around longitude 0°, reaching the 
maximum area at the end of the time series. The formation and 
evolution of individual active regions is mirrored by the beat- 
ings of the amplitude of the light modulation, as can be seen in 
the upper panel of Fig. |2] The minimum amplitude around HJD 
2454285 can be associated with the decay of the active region 
around longitude 0° and the rise of the other one in the second 
active longitude. The second minimum amplitude around HJD 
2454350 is associated with the fading of the spots in both active 
longitudes immediately before the development of the large ac- 
tive region around longitude 0° that dominates the light variation 
until the end of the observations. Therefore, the beating period of 
about 55 days visible in the light curve can be associated with the 
typical lifetime of the largest active regions that form in the two 
active longitudes. A lifetime of about two months is observed 
for the largest and long-lived sunspot groups and compares well 
with those derived from long-term photometry of active stars or 
the models of diffusion of photospheric magnetic fields in stars 
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Fig. 4. Isocontours of the ratio ///max, where / is the spot covering factor and /max = 0.0163 its maximum value, versus time 
and longitude for the ME models with cool spots only. The two dashed vertical lines mark longitudes 0° and 360° beyond which 
the distributions are repeated to easily follow spot migration. The contour levels are separated by A/ - O.l/max with light yellow 
indicating the maximum covering factor and dark blue the minimum. 




Fig. 5. The same as Fig.|4] but for the ME models including spots and faculae with Q - 1.5 and fmax = 0.0203. 



with a low differential rotation amplitude (see, e.g., Ilsik et al] 
l2007l and references therein). 

It is important to note that all individual spots show a remark- 
able retrograde migration after their formation. In other words, 
they have an angular velocity lower than the active longitudes 
by ~ 1.3 percent, as can be derived by the slope of the corre- 
sponding ridges in Fig. |4] This behaviour is reminiscent of that 
of sunspot groups that show an angular velocity about 2 percent 
higher during the first 2-3 days of their lifetime than during the 



later stages of their evolution (e.g., IZappala & Zuccarelloll99l1 : 
Javaraiah & Gokhale 1997). 

The ME models including the effects of faculae show a 
smaller relative range of variation of the spot covering factor and 
a systematic difference in the longitudes of the spot groups ow- 
ing to the peculiar dependence of the facular contrast on the po- 
sition on the stellar disc that makes their photometric effect max- 
imum at some intermediate position between the centre and the 
limb, while spots have a maximum effect closer to the disc cen- 
tre. Nevertheless, the general features of the models with only 
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cool spots are reproduced, although the migration of the second 
active longitude initially at 180° takes place at a slower pace and 
the braking of individual spots is somewhat smaller These re- 
sults indicate that it is difficult to give a quantitative estimate of 
the surface differential rotation or evolutionary effects in the ro- 
tation of starspots in the case of CoRoT-Exo-2a owing to their 
significantly smaller amplitudes with respect to those of the Sun 
and our limited knowledge of the impact of solar-like faculae on 
the observed light modulation. 

The situation is different for the variation of the total spotted 
area, i.e., the integral of the covering factor / over the surface 
of the star It is plotted versus time in Fig. |6] for both the regu- 
larized models with cool spots only and with spots and faculae. 
For comparison, the maximum spotted area observed in the Sun 
is about 0.002 of the surface. In our model, the absolute value of 
the spot coverage depends on the adopted value of the unspot- 
ted magnitude, and the spot and facular contrasts, but its relative 
variations are largely independent of them. In the case of the 
temperature deficit derived for the most active stars, (i.e., adopt- 
ing Cs = 0.21), we find a total spotted area oscillating between 
0.028 and 0.037 of the stellar surface with a variation vs. time 
almost identical to those in Fig. |6l Similar tests show that the 
shape of the modulation of the spotted area is a robust result 
that does not depend on the adopted (uniform) spot contrast or 
the inclusion of faculae in the model. The reported error bars in 
Fig. |6] account only for random errors in the area determination 
and have a semi-amplitude of 3 standard deviations. 

The period of the oscillations of the spotted area obtained 
from the model including only spots is Pqsc = 28.9 ± 4.8 days 
with a confidence level of 99.7 percent (see Fig.|7]l, while for the 
model including also solar-like faculae we find Pose = 29.5 + 4.8 
days with a confidence level of 99.2 percent. Note that the value 
of 28.9 days is close to 10 times the synodic period of the planet 
with respect to the rotation period of 4.522 days of the active 
longitudes. After filtering the variation at the primary frequency, 
we find a secondary periodicity of 86. 5 days, with a conditiona l 
confidence level of only 78.3 percent (iHorne & Baliunaslll986h . 
Its period is three times the principal periodicity and is proba- 
bly associated with the long-term variation of the mean spotted 
area in Fig. |6] although its significance is low given the total 
extension of 142 days of our time series. We investigated the 
impact of temperature variations among different starspots on 
these results by considering two different values of the spot con- 
trast for different longitude intervals. Specifically, we assumed 
a contrast Cs - 0.665 - 6c^ for the longitude interval from 
54° to 234°, roughly encompassing one of the active longitudes 
seen in Figs. |4] and |5] and 0.665 + 6c^ for the remaining lon- 
gitude range, where most of the spots of the other active lon- 
gitude fall. We found that the periodicity at ~ 29 days is still 
detectable for 6c < 0.03 with a confidence level larger than 95 
percent. This corresponds to a relative variation of the contrast 
up to +4.5 percent, i.e., about three times the standard devia- 
tion of the observed sunspot contrast at a wavelength of 672 
nm ( Wesolowski et al. 2008). Note that the relative differences 
among the contrasts of individual sunspots at a given phase of 
the solar 1 1-yr cyc le are gr eater at longer w avelengths (cf., e.g., 
iPenn & Livingston e 2006; P enn & MacDo nald 2007). However, 
solar monochromatic observations at 672 nm are those falling 
closest to the isophotal wavelength of CoRoT observations. 

6. Discussion 

Our spot modelling shows that two active longitudes are present 
on CoRoT-Exo-2a, separated by approximately 180° at the be- 
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Fig. 6. The variation of the total spotted area versus time for the 
ME model with spots only (filled dots; solid line) and the ME 
model including also the effects of faculae with Q - 1.5 (open 
diamonds; dot-dashed line). The area is measured in units of the 
stellar photosphere. 



ginning of the observations. The one showing the largest con- 
centration of spotted area rotates with a period of 4.522 days 
that dominates the light variation as derived from Lomb-Scargle 
periodogram. The other one rotates with a slightly longer period, 
with a relative difference in angular velocity smaller than 1 per- 
cent or even lower if we consider the models including also the 
effects of faculae. The possibility of using active longitudes as 
tracers for the surface differential rotation is supported by the re- 
cent results of Berdyu gina & UsoskinI ( 120031) in the case of the 
Sun, although the propertie s of solar active long itudes are still 
object of debate (see, e.g., lUsoskin et al.ll200"7l for more de- 
tails). The angular velocity of individual spots is smaller than 
that of the active longitudes in CoRoT-Exo-2a, a result that sug- 
gests a dependence of their angular velocity on their evolution- 
ary stage, as observed in the Sun. Those general results are quite 
robust, although a comparison between models computed with- 
out and with faculae shows that solar-like faculae have a sig- 
nificant impact on the determination of the spot longitudes, as 
akeady noticed by iLanza et al.l (|2007|) . Therefore, our results on 
the amplitude of a possible surface differential rotation and the 
spot braking should be taken as upper limits given that they are 
reduced by the facular effects and the true facular area fraction 
Q is not known in the case of CoRoT-Exo-2a. 

It is interesting to note that the braking of sunspots with their 
evolution can be interpreted in terms of a decrease of the depth 
where their flux tubes are anchored, associated with the outward 
decrease of the angular velocity of the Sun in the upper layers of 
the solar convection zone at low la titudes (Schiissler & Rempe| 
I2005t ICorbard & Thompsonll2005h . In other words, our results 
suggest the existence of a subphotospheric shear layer with 
dQ./dr < also in CoRoT-Exo-2a. 

On the other hand, the cyclic variation of the spotted area 
does not depend on the inclusion of the faculae in the model and 
its period can be considered as a robust result. Looking at the 
amplitude of the light curve modulation in Fig. |2] it is not pos- 
sible to see directly the period of 28.9 days because two large 
spots, usually on opposite hemispheres, are modulating the stel- 
lar flux with different amplitudes. However, the periodicity is 
readily apparent when we plot the flux deficit, i.e., the difference 
between the maximum flux and the instantaneous flux, averaged 
over the time along individual rotation periods of 4.52 days (see 
Fig. [8]l. Three cycles are clearly detectable by eye and indicate 
that the oscillation of the spotted area is related to an oscillation 
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Table 2. Parameters of CoRoT-Exo-2a and of the comparison 
stars adopted to check for instrumental effects 



CoRoT ID 


V 


B-V 




P 


U 










(deg) 


(s) 


CoRoT-Exo-2a 


12.568 


0.854 


0.0700 


0.0 


512, 32 


101306449 


12.788 


0.862 


0.0214 


0.225 


512 


101057745 


12.690 


0.938 


0.0359 


0.271 


512, 32 


101297371 


13.310 


0.912 


0.0654 


0.207 


512 


101275308 


12.828 


1.478 


0.1712 


0.111 


512, 32 



of the average light amplitude of about 0.004 mag. We have in- 
vestigated whether such a periodicity may be the result of some 
instrumental effect by performing a Fourier analysis of the sig- 
nals of the brighter stars in the asteroseismic fields, but none has 
been found. Moreover, we have repeated the analysis leading to 
the plot in Fig. [8]for four other stars in the same exofield CCD 
of CoRoT-Exo-2a, after removing long-term trends in their light 
curves with a low-order polynomial, and we found that none of 
them has the same periodicity. Their CoRoT ID, visual magni- 
tude y, color index B-V from Exodat, fraction of photometric 
mask contamination Cm, angular distance p from CoRoT-Exo- 
2a, and exposure time are listed in Table |2] together with the 
same parameters for our tar get. An independen t analysis based 
on the method developed bv iReegen et al.l ( 1200 8) confirms this 
result. Therefore, the cycles of the total spotted area with a pe- 
riod of 28.9 days are likely to be real and not related to some 
instrumental effect. 

The oscillation of the total spotted area with such a short pe- 
riod ma y be regarded as an analogue of the so-called solar Rieger 
cycles. iRieger et al.l d 1 9841) reported a periodicity of about 154 
days in the occurrence of major solar flares close to the max- 
imum of solar cycle 21 and a similar periodicity was found in 
the total sunspot area close to the maxima of cycles 19 and 21 
dOIiver et aI.III998l) . lLoul (l2000l) proposed a theoretical interpre- 
tation based on Rossby-type waves propagating in the subpho- 
tospheric layers of the Sun. Their period is proportional t o the 
rotation period of the star (cf., e.g., Eqs. (13) and (15) in 'LoiJ 
[2OOO). Therefore, the model predicts a period of about 28 days 
in a star rotating about five times faster than the Sun, such as 
CoRoT-Exo-2a. 

An alternative interpretation can be based on the coincidence 
of the period of 28.9 days with 10 times the synodic period of 
the planet with respect to the rotation of the active longitudes. 
This suggests a possible sta r-planet intera ction along the lines of 
the conjecture proposed bv iLanzal (l2008h . Specifically, if some 
subphotospheric dynamo action takes place in the star, the am- 
plified field emerges once th e radial g radient of its intensity 
reaches some threshold value (lAchesonl|T979) . The emergence 
of the flux might be triggered by a small localized perturbation 
of the dynamo or-effect associated with the planet (l Lanzall2b08l) 
when it passes by the active longitudes. An integer number of 
synodic periods is necessary for the planet to trigger again field 
emergence in the active longitudes because the planetary pertur- 
bation can play a role only when the field is already close to 
the threshold thanks to the steady amplification provided by the 
stellar dynamo. 

Unfortunately, the time resolution of our spot maps is only 
3.15 days, i.e., longer than the synodic period, making it impos- 
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Fig. 7. Lomb-Scargle periodogram of the spotted area variation 
for the model including only cool spots (cf. Fig. |6] filled dots). 
The highest peak marked with A corresponds to the periodicity 
of 28.9 days, whereas the lower one, marked with B, corresponds 
to the periodicity of 86.5 days that stands out after filtering the 
modulation corresponding to the principal peak. The horizontal 
dashed line indicates the 99 percent confidence level. 
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Fig. 8. The relative flux deficit versus time, averaged along indi- 
vidual rotation periods of 4.5221 days. 



sible to resolve the effects of individual passages of the planet 
over each active longitude. Making the interval Aff smaller im- 
plies a loss of longitude resolution and the appearance of size- 
able fluctuations of the spotted area when the modelled flux 
modulation covers less than half a stellar rotation. Therefore, we 
must limit ourselves to investigate effects on time scales longer 
than one synodic period. The application of eclipse mapping dur- 
ing successive transits may help to investigate the star-planet in- 
teraction in more detail and will be the subject of a future inves- 
tigation. 

It will be interesting to observe CoRoT-Exo-2a from the 
ground to look for indication of star-planet magn etic interaction 
in the chromosphere, as recently suggested by IShkolnik et al.l 
(2008) for some stars accompanied by hot Jupiters. 

Doppler imaging techniques can in principle be applied 
to CoRoT-Exo-2a because the rotational broadening of photo- 
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spher ic lines as measured by v sin / is ~ 10 km s~^ dBouchv et alj 
l2008 f). Assuming a local line width of 2 - 2.5 km s"' as a con- 
sequence of micro- and macro-turbulence, we expect to have 
4-5 resolution elements across the stellar disk which can pro- 
vide a useful map with some resolution also in the latitudinal 
direction. However, it will be impossible to reach the uniform 
and extended time sampling of CoRoT photometry from the 
ground, even with a coordinated multisite spectroscopic cam- 
paign. Nevertheless, an independent measure of the surface dif- 
ferential rotation with information on the variation of the angu- 
lar velocity v ersus latitude co uld be obtained with the techniques 
developed bv lPetit et all (120021) . 

7. Conclusions 

Our spot modelling of CoRoT-Exo-2a has provided several in- 
teresting results. Solar-like faculae have an area ratio probably 
significantly smaller than in the case of the Sun with Q ~ 1.5 
giving the best fit to the light curves. Spot modelling suggests 
that the differential rotation on CoRoT-Exo-2a is significantly 
smaller than on the Sun. Adopting the relative migration of the 
active longitudes as a tracer for the surface shear, we find an 
amplitude of the relative differential rotation of ~ 0.7 percent 
when only cool spots are considered and even smaller when both 
spots and faculae are included in the spot modelling. Individual 
spots show a rotation rate about 1 percent smaller than that of 
the active longitudes where they form, a phenomenon that can 
be interpreted as due to the presence of a radial negative gradi- 
ent of the angular velocity in the subphotospheric layers and the 
progressive dynamical decoupling of their magnetic flux tubes 
from the deep layers once spots are formed. A quantitative es- 
timate of the amplitude of the spot braking is made difficult by 
our ignorance of the effects of the faculae on the derived spot 
longitudes, specifically by the lack of independent information 
on the Q parameter A robust result, that does not depend on our 
model assumption, is the oscillation of the total spotted area with 
a period of ~ 28.9 days during the five months of CoRoT obser- 
vations. The fact that this period is close to 10 synodic periods 
of the planet with respect to the rotation of the active longitudes 
on the star, suggests some kind of magnetic interaction between 
the star and its hot Jupiter Future studies of the modulation of 
chromospheric indicators, like the flux in the cores of Ca II H & 
K lines, may provide further support for this interaction. We note 
that the detection of a short-term periodicity in a young late-type 
main-sequence star other than the Sun is a new result, since such 
short-term (or Rieger) cycles have been reported only in the RS 
C Vn binary system UX Arietis, consisting of more evolved stars 
dMassi et alj|2005[) . 

The present investigation illustrates the capabilities of 
CoRoT for studying the photospheric activity of late-type stars. 
By applying spot modelling to a large, statistically significant 
sample of targets, we shall be able to understand the evolution of 
their surface magnetic fields on time scales up to a few months 
and use surface brightness inhomogeneities as tracers to study 
differential rotation, opening a new way in the research field of 
solar-stellar connection. 

Acknowledgements. The authors are grateful to an anonymous Referee for a 
careful reading of the manuscript and several valuable comments. This work 
has been partially supported by the Italian Space Agency (ASI) under contract 
ASI/INAF 1/015/07/0, work package 3 170. Active star research and exoplanetary 
studies at INAF-Osservatorio Astrofisico di Catania and Dipartimento di Fisica e 
Astronomia dell'Universita degli Studi di Catania is funded by MIUR {MinLstero 
deU'htruzione, deU'Universitd e della Ricerca), and by Regione Siciliana, whose 
financial support is gratefully acknowledged. This research has made use of the 
ADS-CDS databases, operated at the CDS, Strasbourg, France. 



References 

Acheson, D. J., 1979, Sol. Phys., 62, 23 

Alonso, R., Auvergne, M., Baglin, A., Ollivier, M., Moutou, C, et al. 2008, 
A&A, 482, L21 

Auvergne, M., Boisnai'd, L., Lam-Trong T, et al. 2008, A&A, submitted 

Baglin, A., et al. 2006, 36th COSPAR Scientific Assembly, 36, 3749 

Barnes, J. R, Collier Cameron, A., Donati, J.-E, James, D. J., Marsden, S. C, 

Petit, P 2005, MNRAS, 357, LI 
Berdyugina, S. V. 2005, Living Rev. Sol. Phys., vol. 2, no. 8 
Berdyugina, S. V., Usoskin, L G., 2003, A&A, 405, 1121 

Bouchy, F, Queloz, D., Deleuil, M., Loeillet, B., Hatzes, A. R, et al. 2008, A&A, 
482, L25 

Burrows, A., Budaj, J., Hubeny, I. 2008, ApJ, 678, 1436 

Chapman, G. A., Cookson, A. M., Dobias, J. J. 1994, ApJ, 432, 403 

Corbai-d, T., Thompson, M. J., 2005, Sol. Phys., 205, 21 1 

Donati, J.-E, CoUier Cameron, A., Petit, R 2003, MNRAS, 345, 1187 

Gondoin, P 2008, A&A, 478, 883 

Home, J. H., Baliunas, S. L. 1986, ApJ, 302, 757 

I§ik, E., Schussler, M., Solanki, S. K. 2007, A&A, 464, 1049 

Javaraiah, J., Gokhale, M. H. 1997, A&A, 327, 795 

Kovari, Zs., Strassmeier, K. G., Granzer, T., Weber, M., Olah K., Rice J. B. 2004, 

A&A, 417, 1047 
Lampton, M., Margon, B., Bowyer, S. 1976, ApJ, 208, 177 
Lanza, A. E 2006, MNRAS, 373, 819 
Lanza, A. E 2008, A&A, 487, 1 163 

Lanza, A. F, Catalano, S., Cutispoto, G., Pagano, I., Rodono, M. 1998, A&A, 
332, 541 

Lanza, A. F, Rodono, M., Pagano, I., Barge P, Llebaria, A. 2003, A&A, 403, 
1135 

Lanza, A. E, Rodono, M., Pagano, I. 2004, A&A, 425, 707 

Lanza, A. E, Bonomo, A. S., Rodono, M. 2007, A&A, 464, 741 

Lockwood, G. W., Skifl", B. A., Henry, G. W., Henry, S., Radick, R. R., Baliunas, 

S. L., Donahue, R. A., Soon, W. 2007, ApJS, 171, 260 
Loeb, A., Gaudi, B. S. 2003, ApJ, 588, LI 17 
Lou, Y.-Q. 2000, ApJ, 540, 1 102 

Massi, M., Neidhofer, J., Carpentier, Y., Ros, E. 2005, A&A, 435, LI 

Oliver R., Ballester J.-L., Baudin, F 1998, Nature, 394, 552 

Penn, M. J., Livingstone, W. 2006, ApJ, 649, L45 

Penn, M. J., MacDonald, R. K. D., 2007, ApJ, 662, L123 

Petit, R, Donati, J.-E, Collier Cameron, A. 2002, MNRAS, 334, 374 

Pfahl, E., Arras R, Paxton B. 2008, ApJ, 679, 783 

Pont, F, Gilhland, R. L., Moutou, C, Charbonneau, D., Bouchy, E, et al. 2007, 
A&A, 476, 1347 

Radick, R. R., Lockwood, G. W., Skifl', B. A., Baliunas, S. L. 1998, ApJS, 118, 
239 

Reegen, R, Gruberbauer, M., Schneider L., Weiss, W. W. 2008, A&A, 484, 601 
Rieger, E., Kanbach, G., Reppin, C, Share, G. H., FoiTest, D. J., Chupp, E. L. 

1984, Nature, 312, 623 
Rowe, R. J., Matthews, J. M., Seager S., Kuschnig, R., Guenther D. B., Moffat, 

A. F J., Rucinski, S. M., Sasselov, D., Walker G. A. H., Weiss, W. W., 2006, 

ApJ, 646, 1241 

Rowe, J. E, Matthews, J. M., Seager, S., Miller-Ricci, E., Sasselov, D., Kuschnig, 
R., Guenther D. B., Moftat, A. E J., Rucinski, S. M., Walker, G. A. H., Weiss, 
W. W., 2008, ApJ, in press I arXiv:071 1.41 1 1 1 

Schussler M., Rempel, M. 2005, A&A, 441, 337 

Seager S., Whitney B. A., Sasselov, D., 2000, ApJ, 540, 504 

Shkolnik, E., Bohlender D. A., Walker G. A. H., Collier Cameron, A. 2008, 
ApJ, 676, 628 

Silva, A. V. R. 2003, ApJ, 585, L147 

Solanki, S. K., Unruh, Y. 2004, MNRAS, 348, 307 

Sudarsky, D., Burrows, A., Pinto, P 2000, ApJ, 538, 885 

Usoskin, I. G., Berdyugina, S. V., Moss, D., Sokoloft', D. D. 2007, Adv Space 
Res., 40, 951 

Wesolowski, M. J., Walton, S. R., Chapman, G. A. 2008, Sol. Phys., 248, 141 
Zappala, R. A., Zuccarello, F. 1991, A&A, 242, 480 



